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ABSTRACT 


This thesis presents a computer simulation program designed to predict 
the performance of a laser designator in an operational environment. It 
uses aspects of Fourier optical theory, atmospheric modulation transfer 
function, atmospheric transmittance, and turbulence to analyze the spot size 
and the power on a target. The guidance system requires the _ specific 
power and spot size on a target to track it. So the optimum laser 
designator for special military purpose can be chosen by the output of this 


program for a particular atmospheric condition. 
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I. INTRODUCTION 


A. GENERAL 

In recent years as the naval surface combatants become more 
sophisticated, the use of laser radiation for the detection of targets, and 
weapon control becomes more important for the navy. 

A laser system has several advantages over a microwave system 
because of its very short wavelength. The advantages are higher resolution, 
less power requirement and space needed, interception and jamming are 
difficult, also it can be portable. So, laser guided weapons have been sold 
by the world's weapon manufacturers to all nations capable of supporting 
defense forces. But the laser system has a big disadvantage in poor 
performance in bad weather. 

The performance of military laser systems for imaging, target detection, 
tracking, target designation, warning of missile launch or laser irradiation, 
optical fuzing, or laser weaponry, is very strongly dependent on _ the 
transmission (extinction) and modification of such radiation by atmospheric 
effects. The extinction depends on absorption and scattering by the 
molecular constituents of the atmosphere and by the particles suspended in 
it. The modulation transfer degradation, beam wander, beam _ broadening 
and scintillation effects are caused by atmospheric turbulence. Turbulence 
effects over water as well as over land are a big problem and may be a 


limiting factor in optical systems. 


How to decrease the disadvantage in atmospheric effects and how to 
select the optimum system within the operational environment are more 
important to a user, because of the above factors. So a user (such as the 
Korean navy) must study and research in the area of optical radiation 
propagation through the atmosphere. In support of this effort, special 
attention to computer programming of these effects will be devoted to select 
an optimum laser designator within the operational environment that includes 
the atmospheric transmittance (extinction) and turbulence. This is the 
parameter computed for us by the program. The source of this algorithm 


is based on the program written by Prof. Milne in [Ref. 3]. 


B. BACKGROUND 

For a laser designator, the predicted beam spot size on the target can 
be calculated from the atmospheric modulation transfer function (AMTF) and 
line spread function (LSF) of the laser output pattern and system optics, as 
described below. 

The one-dimensional (single variable) LSF can be generated from the 
two—dimensional point spread function (PSF) by integrating over one of the 
two dimensions of a point source image. This one-dimensional LSF is 
easier to calculate the convolution theorem than the two—dimensional 
function in computer. 

The modulation transfer function (MTF) of the entire system is the 
AMTF times the Fourier transforms of the LSF of the laser output pattern 
and of the optics diffraction (also defocus and aberration broadening). This 


yields the convolution of AMTF, source, and optics. The LSF of the laser 


Output pattern is obtained from its radial distribution by integrating over 
distance y. Then the Fourier transform of the LSF is calculated. 

Taking the inverse Fourier transform of the above yields the LSF for 
the entire system. This LSF can be converted to a radial distribution 
function by application of the ABEL transform. The fractional energy 
falling within a given radius is then calculated by integrating with respect 
to r. 

The algorithm, introduced above, can predict the size of a laser spot 
on a target and the radient energy in the spot for different systems. So, 
we can compare the systems to determine which one is optimum for our 


operational environment. 


ee EO 


The first element in the system we are analyzing is, of course, the 
laser itself. The radial intensity output of the laser is generally Gaussan 


and can be expressed by: 


A = A® exp (-17/207) 


where A; maximum intensity of the laser source in 
watts/steradian 
Tr; angular distance from center in prad 


o; angular divergence of the laser source in prad 


This can be calculated for each laser for approximate results and also can 
be measured directly. One dimensional line spread function for the source 
is obtained from the radial intensity output pattern of the laser by 
integrating to obtain the intensity falling in a vertical slit of width delta x 
(dx) and a distance x from the center of the pattern as shown in Fig 1 
[Ref 1]. A Fourier transform is then used on this LSF to yield the MTF 
of the laser. 

The second element is the optics function. The optics associated with 
the laser is basically near diffraction limited. The Airy function is used to 


calculate the diffraction point spread function for the laser optics. From 


[Ref 6], the Airy function for circular aperture diffraction pattern is 


expressed by the next equation, and illustrated in Fig 2. 


I(@) = I(0) | ee |’ 


where,  —_Ju(u) = first order Bessel function of the first kind 
Ka = hen 
a = radius of the aperture 
E(r) 


(a) Radial distribution of — laser 
output pattern E(r); 


intensity 1; radius 


(b) Integration to produce the line 


spread function intensity vs. x 





is the area of the hatched 


I(x) 
Stri 
* v3 
\ (c) Line spread function intensity 
os vs. X 


(Fig. 1] Change from PSF to LSF 


1/1(0) 


kasin = 


[Fig. 2] Airy function for circular aperture diffraction pattern 


This point spread function is then converted to a line spread function, 
Fourier transformed to get the MTF and then combined by the convolution 
theorem into MTF of the laser plus optics. 

The third element to be calculated is the effects of atmospheric 
turbulance on the laser beam. A_ properly path—position weighted mean 
value of Cy, becomes the central parameter of interest in calculating the 
AMTF due to turbulence. Cp and C; are called respectively the refractive 
index structure constant and the temperature structure constant, and are 


related through the equation below.[Ref 1] 
42 = ( BP 10s oe 


C2 - 2/3 = < (Tan — Tr)? > 


where Try—T ro ; is the temperature difference between two points 
seperated by a distance r 


Ie ; pressure in millibars 


But in the field, the C,y2 can be measured directly by optical methods 
decribed in [Ref 1]. The Cy? is used to predict the MTF of the 
atmosphere over the path between the laser system and the target. the 
AMTF is detailed in subsection 3. 

The AMTF is then multiplied by the transmittance (T = el R) meee 
= extinction coefficient , R = range ) 

Now the total effect of the atmosphere, laser itself and optics can be 
found by the product of the three transform functions.(convolution theorem) 
The inverse transform then will yield the line spread function of the total 
system, the LSF as predicted on the target. [Ref 1] describes the process 
of converting the one dimensional LSF to a circular symmetric point spread 
function (PSF) using the ABEL transform, detailed in subsection 2. 

After obtaining the radial intensity distribution, the size of a laser spot 
on a target can be predicted. Also the power inside a circle of radius r 


can then be obtained by integrating the PSF out to r. 


A. FOURIER TRANSFORM AND CONVOLUTION 
In order to get the total effect of the atmosphere, laser itself and 
optics, some elements of Fourier theory must be considered. The Fourier 


transform represents the one dimensional position variable g(y) expressed in 


the spatial frequency domain u(v). The inverse transform repeats the 
operation in the opposite direction. The forms are commonly represented as 


shown. 


u(v) = Paw exp(—2m ivy) dy = F{g(y)] 


g(y) = [ww exp(2m ivy) dv = F [u(v)] 


The program calculates the value of the integral at a preset number of 
points and yields a discrete Fourier transform. A theorem from Fourier 
theory that allows us to compute the line spread function for one of the 


elements in the system given is the convolution theorem stated as follows. 


If Fig(y)] = G(v), and 
F[h(y)] = H(v) 
then Flg*h] = GH (* represents the convolution) 


B. ABEL TRANSFORM [Ref.1] 

In order to get circularly symmetric two-dimensional PSF from 
one-dimensional LSF the ABEL transform [Griem, 1964; Magnus, 1949] is 
used in this program referred to above in section Il. The one dimensional 
line spread function can be utilized to calculate the corresponding radial 
distribution of intensity in the two—dimensional image, provided that the 
intensity is known to be circularly symmetric, i.e. the intensity E(r) in the 


image depends only on the radial distance from the center of the image. 


This inversion is called the ABEL transform. For such an image, the line 
spread function (LSF) is obtained from the radial distribution by integrating 


the intensity lying within an infinitesmal slit of width dx. 


iz E(r) dy 


I(x) 
R 
=2| E( x? + y? ) dy 


where R is a radius beyond which there is no intensity 


The ABEL transform allows us to reverse this process. It is accomplished 


by carring out the calculation indicated by; 


a 
Ee) at 


where LCS) etegsc Samat 


I’(x) is the x derivative of I(x) 


This integral can be carried out by step wise numerical integration with a 


computer. 


C. MTF OF ATMOSPHERE 
The modulation transfer function is a function of the spatial frequency, 


angular spatial frequency f given in (cycles/radian) or linear spatial 


frequency v given in (cycles/meter). The one-dimensional MTF is directly 
applicable. Images are, of course, really two—dimensional, and much work 
has been done in the field of two—dimensional image processing, both analog 
and digital. | However, for many purposes, the one-dimensional form is 
adequate, and is much more readily handled analytically and in computer 
processing. We already saw that for circularly symmetric images, the 
two—dimensional behavior can be obtained from the one-dimensional. The 
modulation transfer function has the added advantage that for a total 
optical system the system MTF is the product of the MTF's of the 
components. So, the MTF of atmosphere can be expressed independently by 
its own parameters. 

There are two limiting situations used in calculating the AMTF. 
First, the short exposure, or short term MTF is obtained by taking the 
average of successive LSF in which each LSF is shifted such that the center 
of the area of each curve has a common origin. Also it discribes images 
obtained in the limit of time intervals sufficiently short that the turbulence 
can be considered frozen. The average short term MTF is then the average 
of the MTF's corresponding to many such images. This procedure, in 
effect, removes the image wander. Second, the long exposure, or long term 
MTF is obtained by including the image wander, by not performing a shift 
in the origin. Also it results simply from taking an integrated image of a 
sufficiently long time which sees effectively all possible turbulence 
configurations. It is seen to be the result of image spread (the short term 
MTF) plus image wander. So the short term LSF is narrower than the 


long term LSF. 
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From the [Ref 1] AMTF can be summarized by the following 


expression. 


M(f) = exp |-57.64 B Ca? Zo dA 9/3 faq (AR_ps/ ft 


where A = wavelength 
f = angular spatial frequency(cycles/radian) 
D = diameter of optics 
Zo = range 
a = 0 for the L.T. MTF 


= 1/2 for the far-field S.T. MTF, D<<(za)!/2 
= 1 for the near-field $.T. MTF, D>(Zd)!/2 
G6 = 1 for a plane wave 


= 3/8 for a spherical wave 


Di 


Ill. PROGRAM DESCRIPTION 


This program can be run on IBM personal computers or compatibles 
under the DOS operating system (version 2.0 or later), and consists of 
several subroutines, a function and a mainbody. All program modules are 
coded in MS's ANSI 77 Fortran Compiler. The source file is compiled and 
linked to library with command (FL /G2), then the .EXE file is obtained 
and this program is ready to be run. 

The source file is listed in Appendix A, the typical input data file and 
module interconnections are listed in Appendix B. 

some input data and modules are added to the program written by 
Prof. Milne in [Ref. 3] for some parameters to be varied, and for the 


desired output to be taken. 


A. INPUT DATA 
The general input data to this program must be entered interactively 
from the keyboard or from the input file. In either case, these data must 
be determined by the user prior to running this program. The input data 
are described and determined briefly below. 
1. IDATE 
Time and date of the run, 60 characters, it will be placed in the 


header of the output file. 


AZ 


2 INPT 
This input indicates the type of laser source function to be used 
in the program. If 1 is entered, then a gaussian source function is 
calculated. If 2 is entered, then the source function is inputted from the 
input file. If 3 is entered, then combination function of source—optics is 
entered from the input file. 
3. IHEAD(1) 
Information about the laser source, 80 characters, it will be placed 
in the header of the output file. 
4. WAV 
Wavelength of the laser in meters, Format E10.4, this input can 
come from the characteristic of the laser designator determined by the 
weapon manufacturer. 
5. OBRA 
Obscuration ratio, Format F7.5, the ratio of the diameter of the 
inner mirror to outer mirror. This data can come from manufacturer also. 
6. RANGE 
The distance between the laser designator and the target in 
meters, Format F7.0. 
i OPT 
Diameter of the optics in meters, Format F7.3. 
8. SQCN 
Squared refractive index structure constant in meters\-2/ 3) Format 
E10.4. This input will be the path—weighted integrated value in the 
Operational environment and can be measured directly by optical methods 


[Ref. 1]. 
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9. SCALE 
Angular scale in microradians per point, Format F7.0. 
10. EXT 
The extinction coefficient of the operational atmosphere in 
meters!, Format E10.4, this data can come from Lowtran 6 in NPS IBM 
3033AP computer. 
11. INOP 
Input type of the optics diffraction pattern to be used in the 
program. If 1 is entered, then the optics diffraction function is calculated. 
If 2 is entered, then it is entered from the input file. 
12. IHEAD(2) 
Information about the optics, 80 characters, it will be placed in 
the header of the output file. 
13. A 
The maximum intensity of the laser source in watts/steradian, 
Format E10.4, this data can come from manufacturer also. 
14. B 
Sigma(standard deviation) of the laser source (gaussian) in 
microradians, Format E10.4.. This data can be determined to measure the 
angle between the beam center and the line at which the laser's intensity 
falls to e/2 of its peak intensity. This data can come from manufacturer 
also. 
13.7 ee® 
The type of exposure pattern to be taken for output. If 1 is 
entered, then the long—term (exposure) is taken. If 2 is entered, then the 


short—term (exposure) is taken. 
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16. FRAP 
Fraction of total power in watts that is required at the target. 
17. VARY 
Selection of the run. If 1 is entered, then one of input 
parameters will be varied in program. If 2 is entered, then fixed value of 
input parameters will be used in program. 
18, SEL 
selection of the parameter that will be varied in program. 
SQCN, EXT, RANGE, and OPT can be chosen to be_ varied. The 
maximum and minimum of that value should be entered. Input Format is 


same as above 6, 7, 8, and 10. 


B. MODULE DESCRIPTIONS 
1. GETVAR 
a. Module type : Fortran 77 Subroutine 
b. Usage : Call GETVAR(min,max,va) 
Cee Ul GhL Ona To get 20 points between minimum and 
maximum 
d. Interface : 
(1) External references : Float 
(2) Common data : None 
(3) Formal arguments : 
(a) Min (real), minimum of the variable. 
(b) Max (real), maximum of the variable. 
(c) Va (real array), 20 values of the variable between 


minimum and maximum. 
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(4) Input : Min, Max 
(5) Output : Va 
2. AIRY 
a. Module type : Fortran 77 Real function 
b. Usage: Y = AIRY (X) 
c. Function : Calculates the value of the Airy function of x. 
AIRY(x) = 2 * Jy(x)/x, where Ji(x) is the Bessel function of order one. 
d. Interface : 
(1) External references : Float, ABS, MAX] 
(2) Common data: None 
(3) Formal arguments : 
(a) X (real), value of the independent variable of the 
Airy function 
(4) Input : X 
(5) Output : Airy 
3. - LSE 
a. Module type : Fortran 77 Subroutine 
b. Usage : Call LSF(RG3,k) 
c. Function : Converts a point spread function into a line 
spread function. 
d. Interface : 
(1) External references : Float, SQRT, IFIX 
(2) Common data: None 
(3) Formal arguments : 
(a) RG3 (real array), this is the PSF on input and 


the LSF on output. 
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(b) K (integer), number of points in the array RG3. 
(4) Input : RG3, K 
(5) Output : RG3 
4. AMTF 

a. Module type : Fortran 77 Subroutine 

b. Usage: Call AMTF(j) 

Ce Munctiony = Calculates the atmospheric MTF from _ given 
value for long-term and short-term, and then the result is multiplied by 
extinction factor (eHRy 

d. Interface : 

(1) External references : EXP, COS 

(2) Common data: /real/ RG1, RG2. / / WAV, OPT, 
OBRA, SCALE, RANGE, SQCN, EXT, A, B, IDATE, IHEAD 

(3) Formal arguments : 

(a) J (integer), if zero on input, then the long term 

AMTF is calculated; if one on input, then the short term AMTF is 
calculated. 

(4) Input : J, SQCN, RANGE, WAV, EXT, SCALE, 
OPT 

(5) Output : RG2 

oo EEE 

a. Module type : Fortran 77 Subroutine 

b. Usage : Call FFT(RG,ID,KA) 

c. Function : Calculates the Fourier transform or _ inverse 
Fourier transform of a real—value symmetric function using the Fast Fourier 


transform (FFT) algorithm. 
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d. Interface : 
(1) External references : Float, CMPLX, CEXP, CONJG 
(2) Common data: None 
(3) Formal arguments : 

(a) RG (real), on input, RG is an array of functional 
values; on output, it is the forward (or inverse) Fourier transform of these 
functional values. 

(b) ID (integer), specifies whether forward (ID = 1) or 
inverse (ID = 0) Fourier transform is to be calculated. 

(c) KA (integer), the largest non—zero term of forward 
(or inverse) Fourier transform. 

(4) Input : RG, ID 
(5) Output : RG, KA 
6.  BITINV 
a. Module type : Fortran 77 Subroutine 
b. Usage : Call BITINV(K,L,M) 
c. Function : Finds the bit inversion of an M bit number K 
to L. This is used in the Fourier transform subroutine. 
d. Interface : 
(1) External references : None 
(2) Common data : None 
(3) Formal arguments : 

(a) K (integer), the integer whose bit inverse is to be 

found. 


(b) L (integer), the bit inverse of K. 
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(c) M (integer), the number of bits in the integers K 


ama: 1. 
(4) Input : K, M 
(5) Output : L 
7. ABEL 


a. Module type : Fortran 77 Subroutine 

b. Usage : Call ABEL(K) 

c. Function : Converts a line spread function (x) to a point 
spread function (r). 

d. Interface : 

(1) External references : Float, SQRT, ABS 
(2) Common data: /real/ RG1, RG2 
(3) Formal arguments : 
(a) K (integer), number of points in array RG2 to be 
processed. 
(4) Input : RG2, K 
io meOutput a RG? 
8. INT 

a. Module type : Fortran 77 Subroutine 

b. Usage: Call INT(K,TPOW) 

c. Function : Integrates the point spread function to find the 
total power as a function of radius reached around the target and 
normalizes it. 

qd. "Interface : 


(1) External references : Float 
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(2) Common data: /real/ RGl, RG2. / / WAV, OPT, 
OBRA, SCALE, RANGE, SQCN, EXT, A, B, IDATE, IHEAD 
(3) Formal arguments : 

(a) K (integer), number of points in array RG2. 

(b) TPOW (real), total integrated power reached 
around the target. 

(4) Input : K, RG2, SCALE, RANGE 
(5) Output : TPOW 
9. OUTIN 
a. Module type : Fortran 77 Subroutine 
b. Usage : Call OUTIN(BUF,J,K,L) 
c. Function : Control the input and output files to read and 
write the data. 
d. Interface : 
(1) External references : None 
(2) Common data : /OI/ LU, LI, IH, IM, IS, IHS, IY, 
IMO, ID 
(3) Formal arguments : 

(a) BUF (array), array data which are to be read or 
written. 

(b) J (integer), if —l1 on input then call the OPEN 
subroutine; if 0 on input then call the READF subroutine; if 1 on input 
then call the WRITF subroutine. 

(c) K = (integer), the count number of variable to be 


varied, 


(d) L (integer), parameter specifying which title is to 


be printed on output file or which data is to be read from input file. 


10. 


IMO, ID 


ae 


IMO, ID 


data file. 


(4) Input : BUF, J, K, L 
(5) Output : BUF 
OPEN 
a. Module type : Fortran 77 Subroutine 
b. Usage: Call OPEN 
c. Function : Open the input and output files. 
d. Interface : 
(1) External references : None 


(2) Common data : /OI/ LU, LI, IH, IM, IS, IHS, IY, 


(3) Formal arguments : None 
(4) Input : None 
(5) Output : None 
READF 
a. Module type : Fortran 77 Subroutine 
b. Usage : Call READF(IER,BUF,L) 
G: Function : Read the data from input file. 
d. Interface : 
(1) External references : None 


(2) Common data : /OI/ LU, LI, IH, IM, IS, IHS, IY, 


(3) Formal arguments : 


(a) IER (integer), error number of reading the input 


Zi 


(b) BUF (array), array data which are to be read 
from input data file. 
(c) L_ (integer), if 1 on input then read the 


information; if 2 on input then read the source function or source—optics 


function. 
(4) Input: L 
(5) Output : IER, BUF 
12. WRITE 
a. Module type : Fortran 77 Subroutine 
b. Usage: Call WRITF(IER,BUF,L) 
c. Function : Write the data in output file. 
d. Interface : 
(1) External references : None 
(2) Common data : /Ol/ LU, LI, TH, IM, 15, aise 
IMO, ID 
(3) Formal arguments : 
(a) JER (integer), error number of writing the output 
data file. 


(b) BUF (array), array data which are to be written 
on output data file. 
(c) L (integer), parameter specifying which title is to 
be printed on output file. 
(4) Input : BUF, L 
(5) Output : IER 
13. RESUL 


a. Module type : Fortran 77 Subroutine 


Ze 


b. Usage: Call RESUL(INU,TPOW,RADIU,POWER,FRAP) 
c. Function : Calculates the radii on the target and average 
usable power. 
dee lntertaces: 
(1) External references : Float 
(2) Common data: /real/ RGl, RG2. / / WAV, OPT, 
OBRA, SCALE, RANGE, SQCN, EXT, A, B, IDATE, IHEAD 
(3) Formal arguments : 
(a) INU (integer), the count number of variable to be 
varied. 
(b) TPOW (real), total integrated power reached 
around the target. 
(c) RADIU (real array), radii on the target which 
depends on the variable VA. 
(d) POWER (real array), average usable power which 
depends on the variable VA. 
(e) FRAP (real), function of total power that is 
required at the target. 
(4) Input : RG2, FRAP, RANGE, SCALE, TPOW 
(5) Output : POWER, RADIU 


14. WRITO 
a. Module type : Fortran 77 Subroutine 
b. Usage: Call WRITO(VA,RADIU,POWER,ABC) 
c. Function : Write the data in OUT2. file. 


d. Interface : 
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(1) External references : None 
(2) Common data: /OJ/ SEL, VERY 
(3) Formal arguments : 
(a) VA (real array), 20 values of the variable between 
minimum and maximum. 
(b) RADIU (real array), radii on the target which 
depends on the variable VA. 
(c) POWER (real array), average usable power which 


depends on the variable VA. 
(d) ABC (array), information that includes the 


common data / / WAV, OPT, OBRA, SCALE, RANGE, SQCN, EXT, A, 


B, IDATE, IHEAD. 
(4) Input : VA, RADIU, POWER, ABC,SEL 


(5) Output : None 
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IV. USE OF PROGRAM 


The purpose of this simulation is to use the program to predict the 
atmospheric effects on a laser designator, and to select the optimum laser 
designator within the operational environment. This simulation program 
includes only one leg from a laser designator to a target, but the total 
infrared system including laser designator must be described before simulating 
with this program. 

The total infrared system including laser designator with definite 
military mission can be represented by the following block diagram shown in 
Fig. 3. The specific function of the blocks are discussed in [Ref. 7]. 

This total system involves double pass operation from a_ transmitter 
(laser designator) to a target and back to the receiver. The received signal 
may be used to determine the position, range, or velocity of the target, 
then the guidance system can be directed at the _ target. A target 
designator system is the first leg of this total system in which the receiver 
is physically removed from the transmitter. The most obviously important 
measure of system performance is the maximum range of operation, or the 
ratio of received power to transmitted power. The range equation for this 


double pass is stated by following expression in [Ref. 7]. 


|e ie 
— us d m 
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where P, = output power of the transmitter (laser designator) 


in watts. 
P, = __ received power of the detector in watts. 
At = aperture of the transmitter in meters. 
Ag = __ aperture of the detector in meters. 
m= overall path efficiency (= Ta x Tm =x To, Ta is 


atmospheric transmittance from designator to 
target, Im is atmospheric transmittance from target 


to detector, T> is optical transmittance). 


Ra = _ distance of designator from target. 
Rn = _ distance of detector from target. 
o¢ = optical reflection efficiency of the target (the 


product of spot area on target and_ reflection 
efficiency of target surface and cos6). 
This equation based on some assumptions, the designator (transmitter) and 
detector are on different path, so the angle between designator and detector 
on target is 6, the target acts as a diffuse radiator, and the spot size on 
the target is smaller than the target size. 

From the assumptions, the reflected flux per unit solid angle for a 
perfectly diffuse reflector is proportional to the cosine of the angle between 
direction of interest and the normal to the surface, a relationship known as 
Lambert's cosine law.[{Ref. 8] 

For the guidance system shown in lower part of Fig. 3, there exists a 
minimum IR power input from the target necessary for detection. This 
must be in the form of a target irradiance at the optical receiver in 


guidance system. The minimum detectable spectral intensity from the target 
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depends on the distance between the target and the optical receiver, and 
minimum spectral power required in the receiver that includes the optical 
transmission factor and spectral responsivity of detector. So the range 
equation for guidance system can be constructed to get the minimum 


detectable spectral intensity from the target by: 


It(A) = P(A) - Ro 
(Np ee ASUS 
where Irt(A) spectral intensity reflected from target 
(Watts/sr) 
Pr{A) minimum spectral power required in receiver 
(Watts) 
Ao ; aperture of the detector (m2) 
u() ; spectral extinction coefficient (m-!) 
R range between target and guidance system 


Finally the required spectral intensity on target from the laser 


designator can be achieved considering the optical reflection efficiency of 


target. 
Tit(A) = erie. = 
where lit(A) required spectral intensity on target from the 


laser designator (Watts/sr). 


0 optical reflection efficiency of the target. 
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In order to track a target, the tracker or the guidance system must be 
directed at the target within a certain angle, and the signal to noise ratio 
of the target signal must be large enough that the tracker can detect the 
target. Thus the prime military requirement for a laser designator is 
optical power on the selected target. The required spectral intensity on 
target must be achieved in a variety of operating conditions. 

And now, the first leg of this total system shown in upper part of 
Fig. 3 can be treated to select the optimum designator in operational 


environment, and was decribed in Chapter II. 
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A. METHOD OF SELECTION 

The typical values of laser designators and guidance systems are 
classified and are not allowed to public use. For the purpose of this 
simulation, the assumed data is to be used or only the method for test is 
to be explained. 

Also the one sigma (lo) point (40%) of total integrated power is to 
be used to predict the spot size on target place. The value of sigma 
represents the inflection points of the Gaussian laser power output (the 
points at which the curve changes from turning downward to turning 
upward), and 40% of the output power is included within these points. 

The averaged usable power on target is to be used to predict the 
available power for tracking and guidance system. This power is more 
useful considering pointing accuracy and tracking stability. 

Four parameters, Cy2, Extinction coefficient, Diameter of optics, Range 
are available to be varied in this simulation program. 

Now some examples of using this method for this simulation program 
are to be introduced. First, how is the spot size changed when Cp? is 
varied. Second, how is the usable power changed when extinction coefficient 
is varied. Third, how is the spot size changed when range is varied. 
Fourth, how can the system parameter be selected when specific operational 
environment is given (single value of Cy2). Finally, how can the optimum 
system be selected when special operational environment is given (various 
values of Cy2). 

1. TURBULENCE EFFECTS (RUN 1) 

a. Input data 
(1) Wavelength ; 1.06 ym 
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target 


Obscuration ratio ; 0.00 
Diameter of optics ; 0.95x10°! m 
Scale ; 1.0 prad/point 


Max. intensity ; 500 watts/steradian 


Cr? + 8.885x10-13 — 2.258x10-13 m-2/3 
Range ; 10 Km 


Extinction coefficient ; 0.5x10°4 m7! 


(2) 

(3) 

(4) 

(5) 

(6) Sigma ; 1.5 prad. 
(7) 

(8) 

(9) 

( 

( 


b. Output ; see Appendix C, A. 

EXTINCTION EFFECTS (RUN 2) 

a. Input data ; same as RUN 1 except (7) and (9) 
(7) Cph2 ; 1.0x10-14 m-2/3 
(9) Extinction coefficient ; 0.5x10°5 — 0.5x10°3 m7! 

b. Output ; see Appendix C, B. 

RANGE EFFECTS (RUN 3) 

a. Input data ; same as RUN 1 except (7) and (8) 
(7) Cp2 ; 1.0x10714 m-2/3 
(8) Range ; 5 — 15 Km 

b. Output ; see Appendix C, C. 

SELECTION OF DIAMETER (RUN 4) 

a. Input data ; same as RUN 1 except (3) and (7) 
(3) Diameter of optics ; 0.25x10°! — 1.25x10°! m 
(7) Cp2 ; 1.0«10714 m-2/3 
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b. Output ; see Appendix C, D. 


COMPARE TWO SYSTEMS 
Here two individual system parameters like (1) through (6) in 
But the environment parameters like 


D. 


RUN 1 are to be used for each run. 
(7) through (10) in RUN 1 are to be taken in the same values, because 


these two systems must be compared in same atmospheric environment. 


Also real spectral intensity required on target [Iit(A)] instead of 


one sigma point must be used to determine which one is more optimum for 


the purpose including tracking or guidance system. 
As a result the output will show spot sizes and fractional powers 


inside a circle for each run, then the user must make a decision which 


system is the best for their purpose and operational environment. 


B. RESULTS & CONCLUSION 
This program can be useful in predicting the performance of laser 


It was used to show how the spot 


designator in operational environment. 
size and the averaged usable power are changed when some parameters are 


varied as shown in above section A. 
The C,r2 values were taken for 


Cn2 was. varied. 


In Run 1, the 
As a result the radius of spot size was 


area from Ref. 3. 


Monterey bay 
0.24m to 1.36m as well as the averaged usable power was 


changed from 


997 to .2726x10"!watts/sr. 
The radius of spot size 


was changed from 


changed from 
2, extinction coefficient was varied. 


the averaged usable power 
From the results of Run 1 and Run 2, 


In Run 
was not changed but 
.1411x10*! to .9998x10°2 watts/sr. 
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it is shown that the performance of a laser designator is affected 
significantly by atmospheric effects. 

In Run 3, range was varied. As a result the radius of spot size was 
changed from .09 to .465m as well as the averaged usable power was 
changed from 2.234 to .4551 watts/sr. From the results, the maximum and 
minimum range can be predicted in specific environment. 

It was used to select the optimum diameter of optics as shown in 
Run 4. The variation was not too big, and over .83cm of diameter the 
radius of spot size and averaged usable power was not changed. From the 
results, the optimum diameter of optics can be predicted. 

This program can also be used to select the optimum system in special 
operational environment and purpose like the Pacific Ocean, the Atlantic 
Ocean, and air—to—air, air—to—surface, surface—to—surface systems. 

In the military use, designators are used to mark or designate a 
target, thereby turning the guidance system into a beacon of reflected 
energy from target. Given this use, the required power on the target must 
be achieved over a variety of operating environments. Also laser system 
user must look at such factors as target size and target cross section, 
pointing accuracy and tracking stability requirement, as well as atmospheric 
conditions like turbulence and attenuation. Designators will have energy 
amplitude stability, sidelobe, and radiation outside the main beam 
requirements to avoid loss of lock—on or misdirection of seekers sensing the 
reflected energy. Also no visible light can be emitted and the beam width 


must be reduced to minimize countermeasure effectiveness. 
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APPENDIX A 


PROGRAM LISTING 


Meee C COCOA i a 
C 


C 

C THIS PROGRAM PREDICTS THE PERFORMANCE 

C OF THE LASER DESIGNATOR WITHIN 

C OPERATIONAL ENVIRONMENT TO SELECT THE OPTIMUM SYSTEM. 
C 
C 
C 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


INTEGER*2 IDATE, IHEAD, INFO,IH,IM,IS,IHS,IY,IMO,1ID,LI,LU,SEL, VARY 
COMMON /REAL/RG1 (512) ,RG2(512) 


COMMON//WAV, OPT, OBRA, SCALE, RANGE, SQCN ,EXT, A,B, IDATE(30) , IHEAD(80) 


COMMON /O1/LU, LI, IH, IM, 1S, IHS, IY, IMO, ID 
COMMON /0)/SEL VARY 
DIMENSION INFO(128). ,ABC(64) ,VA(20) ,RADIU(20) , POWER (20) 
EQUIVALENCE (WAV, INFO(1)) , (INFO(1) ,ABC(1)) 
( DATA INPUT SECTION 
LI=7 
LU=8 
CALL GETDAT(IY, IMO, ID 
CALL OUTIN(ABC,-1,1,1 
CALL OUTIN(ABC,0,1, 1) 
WRITE(*, 1000) (IDATE(1) , I=1,30) 


1000, FORMAT(' ENTER TIME AND DATE OF THE RUN - 60 CHARACTERS " /TXSSUR2 


EAD(*, 1030) (IDATE(T) ,1=1,30) 
10 WRITE(*,1010) 


C 
C 
C 
C 
C 
C 
C 


1010 FORMA Pt 1 IF SOURCE IS GAUSSIAN TO BE CALCULATED'/' ENTER 


$2 IF TH 
$FUNCTION IS TO BE INPUTTED _‘'') 
READ(* .*) *)INPT 
WRIT iy 1020) (THEAD (I I=1,40) 
1020 FORMAT(/' R SOURCE INFORMATION - 80 CHARACTERS'/1X,40A2, /) 
READ (*, 1030 (IHEAD (I), I=1,40) 
1030 FORMAT (4042 
WRITE(*,1040)WAV 
1040 dite EN ER WAVELENGTH IN METERS! ,2X,E10.4, /) 
READ 
WRITE. 105 1050 OBRA 
1050 FORMAT (/' ENTER OBSCURATION RATIO',2X,F7.5,/) 
READ (*, *)OBRA 
WRITE (*, 1060) RANGE 
1060 FORMAT(/! ENTER RANGE IN METERS' ,2X,F7.0, /) 
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OURCE IS TO BE INPUTTED'/' ENTER 3 IF THE SOURCE- OPTICS 


ead | *) RANGE 


WRIT (, ,1070) OPT 

1070 FORMAT (/' ENTER DIAMETER OF OPTICS IN METERS! ,2X,F7.3,/) 
READ (*,*) OPT 
mere 1080) SQCN 

1080 FORMA iat CN SQUARED IN METERS**(-2/3)',2X,E10.4, /) 

EAD 

weITEG 1090 SCALE 

1090 FORMA (/' ENTER SCALE IN MICRORADIANS PER POINT! ,2X,F7.2,/) 
READ (*,*) )SCALE, 
WRITE(*, 1100 

1100 FORMAT () EN rR EXTINCTION COEFFICIENT IN (1/METERS)',2X,E10.4, /) 
READ (* *)EXT 


1F (I PT.£Q.3) GO 70 11 
WRITE(*, 1015 
1015 FORMA ( ' ENTER 1 IF OPTICS DIFFRACTION PATTERN IS TO BE CALCULATE 
$D'/! ENTER 2 IF OPTICS DIFFRACTION IS TO BE INPUTTED") 
READ (* ,*) INOP 
WRITE(*, 1110) (THEAD(T) I=41, 80) 
1110 FORMAT(/' ENTER OPTICS INFORMATION - 80 CHARACTERS' /2X,40A2, /) 
READ (*, 1030) ( (THEAD (I), 1=41,80) 
11 IF (INPT-2) 12,20,65 
C CALCULATION OF SOURCE PATTERN IF NOT INPUTTED 
12 WRITE(*,1120)A 
1120 age ( ENTER MAXIMUM OF SOURCE IN WATTS/STERADIAN ',2X,E10.4,/) 
READ (* 
WRITE bee B 
1130 FORMA (/' ENTER SIGMA OF GAUSSIAN SOURCE IN MICRORADIANS' ,2X,E10.4 


*} gh | 
WRITE(*, 1140 
1140 FORMAT(/' ENTER 1 IF LONG- TERM (EXPOSURE) IS TO BE TAKEN FOR OUTPUT 
s" ' ENTER 2 IF SHORT-TERM(EXPOSURE) IS TO BE TAKEN FOR OUTPUT') 
READ(*,*) EXPO 
WRITE(*, 1150) 
1150 FORMAT (/ ENTER THE FRACTION OF TOTAL POWER THAT IS REQUIRED AT TH 
SE TARGET 
READ(*,*) FRAP 
WRITE(*, 1160 
1160 ri ENTER 1 IF ANY VARIABLE IS 10 BE VARIED'/' ENTER 2 IF AN 
$Y ae NOT TO BE VARIED ') 
READ (* ,*) VARY 
eee} 110,110,100 
110 WRITE(*, 1170 
1170 PORMAT(/" 4 | ARTABLES CAN BE VARIED IN THIS PROGRAM /) SELECT THE 
HAT 


SNUMBER IS TO BE VARIED'//' 1 - SQCN'/' 2 - EXT'/' 3 - RANGE'/ 
a) ae ee eTER) ) 
READ (*,*) SEL 


GOTO (120,130,140,150) SEL 
120 WRITE(*, 1180) 
1180 FORMAT(/' ENTER THE MIN-VALUE OF SQCN') 
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1190 


200 
130 
1200 


1210 


300 
140 
1220 


1230 


400 
150 
1240 


1250 


900 
100 


1260 


14 
15 


READ (*, *) SQCNL 
WRITE(*, 1190 
FORMAT(/' ENTER THE MAX-VALUE OF SQCN') 
READ(*, ys CNH 
CALL GETVAR(SQCNL,SQCNH,VA) 
INU=1 


U= 
SQON-VA(INU) 
TO 100 
WRITE 
FORMAT (/! UNTER THE MIN- VALUE OF EXT’) 
READ (*, *) EXTL 
WRITE(*, 1210 
min ENTER THE MAX-VALUE OF EXT') 
READ (*, *) EXTH 
CALL CORT ARCEXTL, EXTH, VA) 
INU=1 
EXT=VA (INU) 
GO TO 100 
WRITE 
PORMAT( 7! "ENTER THE MIN- VALUE OF RANGE") 
Tg 
WRITE (*, 1230 
FORMAT(/' ENTER THE MAX- VALUE OF RANGE" ) 
READ (*,*) RANGH 
CALL BTV AR(RANGL, RANGH , VA) 
INU=1 
RANGE=VA (INU) 
TO 100 
WRITE(*,1240 
FORMAT (/" ENTER THE MIN- VALUE OF OPT(DIAMETER) ') 
READ (* ;")OPTL 
WRITE ( 
FORMAL (]" "oxo THE MAX- VALUE OF OPT(DIAMETER) ') 


CALL GbTVAR(OPTL, OPTH,VA) 
INU=1 

OPT=VA (INU) 
IF (VARY.EQ.2) INU=20 

WRITE(*, ,1260) INU,VA(INU) 

FORMAT(/' CONTINUE ,2X,I2,2X,E10.4) 
DO 15 [=1,512 

F=SCALE*FLOAT (I- 1) 

ARG=F*F / a _*B*B) 


IF (ARG. G at TO 14 
RGY (1) APEXPL- ARG) 

K=I 

GO TO 15 

ROL (I)= 

CONTINUE 

CALL OUTIN(RG1,1,INU,2) 
GO TO 25 
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Ca'ea 


INPUT OF SOURCE PATTERN IF NOT CALCULATED 
20 CALL OUTIN(RG1,0,20,2) 
DO 22 I=1,512 
a wii 
IF(RG1(K) .NE.O.) GO TO 25 
22 CONTINUE 
SOURCE PSF IS CONVERTED TO LSF 
25 CALL LSF(RG1 K) 
CALL OUTIN(RG1,1, INU, 10 
FOURIER TRANSFORM OF SOURCE LSF 
CALL FPT (RG 1,K) 
CALL OUTIN(RG1,1,INU, 18) 
IF (INOP- 1)30,30,40 
CALCULATION OF OPTICS DIFFRACTION PATTERN IF NOT INPUTTED 
30 OBS=0BRA*OBRA 
OB=(1.- 0BS)**2 
Q=3.14159E- 6*SCALE*OPT/WAV 
DO 34 I=1,512 
R=Q*FLOAT(I- 1) 
IF(R.GT. 60. ne 
teri) =((AIRY(R 


00 TO 34 
32 RO2(1)=0. 
34 CONTINUE 
CALL OUTIN(RG2,1, INU, 26) 
GO TO 46 
INPUT OF OPTICS DIFFRACTION IF NOT CALCULATED 
40 CALL OUTIN(RG2,0,20,26) 
DO 44 I=1,512 
K=512-1 
1F(RG2(K). .NE.O.) GO TO 46 
44 CONT 
CONVERSION OF OPTICS DIFFRACTION PATTERN (PSF) TO OPTICS LINE SPREAD 
LSF 
46 CALL LSF(RG2,K) 
CALL OUTIN(RG2,1, INU. 34) 
FOURIER TRANSFORM OF OPTICS LSF 
CALL FFT(RG2,1,K 
NORMALIZATION OF THE FOURIER TRANSFORM OF OPTICS LSF 
DIV=RG2(1) 
DO 57 I=1,512 
57 RG2(1)=RG2(I) /DIV 
CALL OUTIN(RG2,1,INU,42) 
PRODUCT OF FOURIER TRANSFORM OF SOURCE LSF WITH FOURIER TRANSFORM OF 
LSF 
DO 64 I=1,512 
64 RG1(1)=RGi (1) *RG2(I) 
CALL OUTIN(RG1,1,INU,50) 
GO TO 66 
INPUT OF FOURIER TRANSFORM OF SOURCE LSF * FOURIER TRANSFORM OF OPTIC 


)- OBS*AIRY (R*OBRA))**2) /OB 
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65 CALL OUTIN(RGL. ,0,20,50) 
C CALCULATION OF ATMOSPHERIC MODULATION TRANSFER FUNCTION (LONG TERM) 
66 CALL AMTE (0) 
CALL OUTI G2.1, INU. 58) 
C PRODUCT OF AMTF WITH FFT OF SOURCE LSF AND FFT OF OPTICS LSF (LONG T 
DO 68 I=1 it 
68 RG2(I)=RG1(I er ” 
CALL OUTIN(RG2,1,INU,66) 
C INVERSE FOURIER TRANSFORM bE IMAGE LSF (LONG TERM) 
CALL PPT (RG2, 0,K) 
CALL OUTIN(RG2,1,INU,74) 
C CONVERSION OF IMAGE LSF TO PSF (LONG TERM) 
CALL ABEL K) 
CALL OUTI G2,1, INU,82) 
CALCULATION OF FRACTION OF POWER TO TARGET INSIDE A CIRCLE OF RADIUS 


ap eee, 


IF (EXPO. EQ. 2) GOTO 7 
CALL RESUL (INU, POW. RADIU, POWER , FRAP) 
WRITE(*,1270) TPOW,RADIU(INU) , FRAP, POWER(INU) 
1270 FORMAT (/' TOTAL P=',E10.4,2X, ‘RADIUS="' ,E10.4,2X, 'FRAC=',F5.3,2X,' 
SAVE P=", E10.4) 
72 CALL OUTIN(RG2, 1, INU, 90) 
C CALCULATION OF ATMOSPHERIC MODULATION TRANSFER FUNCTION (SHORT TERM) 
CALL raat 
CALL OUTIN(RG2,1, INU, 98) 
C PRODUCT OF ni WITH FFT OF SOURCE LSF AND FFT OF OPTICS LSF (SHORT 
DO 74 I=1,512 
74 RG2(1)=RG1 (1) *RG2(I 
CALL OUTIN(RG2,1, INU, 106) 
C INVERSE FOURIER TRANSFORM OF IMAGE LSF (SHORT TERM) 
CALL FFT(RG2,0,K) 
CALL OUTIN(RG2,1, INU, 114) 
C CONVERSION OF IMAGE LSF TO PSF (SHORT TERM) 
CALL ABEL( (") 
CALL OUTIN(RG2, 1, INU, 1122) 
CALCULATION of FRACTION OF POWER TO TARGET INSIDE A CIRCLE OF RADIUS 
(SHORT TERM) 
CALL INT K TPOW) 
IF (EXPO. EQ. 1) GOTO 76 
CALL RESUL(INU,TPOW, RADIU, POWER, FRAP) 
WRITE(*,1270) TPOW,RADIU(INU) , FRAP, POWER (INU) 
76 CALL OULIN(RG2,1, INU, 130) 
INU=INU+1 
IP(SEL.BQ)1 AND. INU.LT.21) THEN 
0TO 200 
ELSEIF (SEL.EQ.2.AND.INU.LT.21) THEN 
GOTO 300 
ELSEIF (SEL.£Q.3.AND.INU.LT.21) THEN 
GOTO 400 
ELSEIF (SEL.EQ.4.AND.INU.LT.21) THEN 


ag | ae. 
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80 


—~J OT 


20 


10 


30 


40 


20 


60 


fo | 


10 


100 
80 


GOTO 500 
ELSE 


GOTO 80 
ENDIF 
CALL WRITO(VA,RADIU, POWER, ABC) 
CALL OUTIN(ABC, 2,20, 1) 
STOP 
END 


SUBROUTINE RESUL(INU, TPOW, RADIU, POWER, FRAP) 
DIMENSION ae POWER (20 
COMMON /REAL/RG1 (512) ,RG2(512 
COMMON//WAV, OPT, OBRA, SCALE, RANGE, SQCN, EXT, A,B, IDATE(30) , IHEAD(80) 
DO 5 I=1,512 
P_RG2(T) /FRAP 
1F(P. E.1.0) GOTO 7 
CONTINUE 
RADIU( INU) =RANGE*FLOAT (1) *SCALE*1 . E- 6 
POWER (INU) =(TPOW*FRAP) /(3.14159*RADIU(INU) **2) 
ETUR 


END 


SUBROUTINE WRITO(VA,RADIU, POWER, ABC) 
INTEGER*2 SEL, VARY 
DIMENSION VA(20) =RADIU(20), POWER (20) , ABC(64) 
COMNON/0J /SEL.,V 
IF (VARY .EQ 2) OTD 100 
DO 10 I=1 
WRITE(9, 30) 'VA(L), RADIU(I) , POWER(I) 
FORMAT (2X,£10.4,2X,F10.3,2X,E10.4) 
CONTINUE 
WRITE (9,30) 
FORMAT(2X,' VARIABLE',7X, RADIUS! ,3X,'AVERAGE USABLE POWER' ) 
GOTO (2 3, 4,5) SEL 
WRITE(9, 40) 
FORMAT(/,2X,'VARIABLE IS SQCN',/) 
GO TO 200 
WRITE(9 50 
FORMA (. XK. 'VARIABLE IS EXT’, /) 
GO TO 2 
MEITE 360 
PORMAT(/. (J, OK. 'VARIABLE IS RANGE', /) 
T 
WRITE(9, 70) 
FORMAT (/, 2X, VARIABLE IS OPT(DIAMETER) ', /) 
TO 2 
WRITE(9,80) RADIU(20) ,POWER(20) 
FORMAT (2X,F10.3,2X,E10.4, /) 
WRITE(9, 90) 
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90 
200 
300 


10 


FORMAT (2X, "RADIUS! ,6X,'AVERAGE USABLE POWER") 

WRITE(9, 1300) (ABC 1) , I=1,64) 

FORMAT (' WAVELENGTH = ' £12.5,' DIAMETER OF OPTICS = 
$E12.5,/, OBSCURATION RATIO - | 812.5," SCALE - 
§ £12.5,/, 'RANGE = ' £12.5,' CN SQUARED : 
$i £12.5,/, EXTINCTION COEEFICIENT = ' £12.5," AMPLITUDE = 
§ ' £12.5,/, "SIGMA  E12.5,/,2X, 1544, /,2X,20A4, / 
§ 2X, 20A4) 

RETURN 

END 


SUBROUTINE GETVAR (NIN, MAX ,VA) 
DIMENSION VA(20 

REAL MIN,MAX,K, 

K=MAX- MIN 

L=K/19. 

VA(1) =MIN 

DO 10 I 

VA(T) = VA (1) sLFLOAT(I- i) 
CONTINUE 

VA(20 MAX 

RETUR 

END 


SUBROUTINE ABEL (K) 


C CONVERTS LINE SPREAD FUNCTION (X) TO POINT SPREAD FUNCTION (R) 


rr 


GO bo 


oa 


ee RG2 (512) 
G=RG2(1) 

F=G 

ee .2* (RG2(1)- RG2(2) ) 

We eal Se 

H=RG2(I+1) 

RG2(I+1)=.4*G+.2*H- .6*RG2(1+2) 
G=H 

il G).GT.1E-6*F) K=I 
CONTIN 


DO 3 I=1,K 
RG2(1) =R62(T) /(2.*SQRT\ (FLOAT(I)- .9)**2- (FLOAT(T- 1))**2)) 


DOP2 J=IhkK 

ra =RG2(I)+RG2(J) /SQRT( (FLOAT (J)- .9)**2- (FLOAT (I- 1) )**2) 
RG2(1)=RG2(1) /3.14159 

ae Wi 

K=K+1 

DO 4 [=k,512 

RG2 a 

RET 

END 
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FUNCTION AIRY (X 


C ATRY(X) : sth) WHERE J1(X) IS THE BESSEL FUNCTION OF ORDER ONE 
TO 8 


wh —_ 


i CO 


Or 


—~I Oo) 


oe @) 


IF(X.EQ.0.) G0 
z 


ot) GO TO 1 
J=20.4+10. *X- X*X/3 
GO TO 2 

J=90.+X/2. 
IF(X.LT.5.) GO TO 3 
W=1.4*X+60./X 

GO TO 4 

W=X+6. 

T=MAX1(W, (3.+X/4.)) 
DO 6 M=I,J,3 

E- 28 


i= mnamore 


= ee 


nr 
Won Ww ul oe 
a™ Cat <t bS 


M-K-2).EQ.0) B=U 


an 
a 
TT 
poeh 
q 
— 


S=S+V*FF 

U=2. *V/X-T 

S=S+U 

B=B/S 

ee KEE 6)) 17.0.) GO TO 7 


CONTINUE 
ATRY=2.*B/X 
RETURN 
ATRY=1. 
RETURN 

END 


SUBROUTINE AMTF (J) 


C THIS SUBROUTINE CALCULATES THE ATMOSPHERIC MODULATION TRANSFER FUNCTI 
C FOR LONG TERM (J=0) AND SHORT TERM (J=1) AND STORES THE RESULT IN RG2 
C THE RESULT IS THEN MULTIPLIED BY EXTINCTION FACTOR EXP(-EXT*RANGE) . 


COMMON /REAL/RG1 (512) ,RG2(512) 
COMMON/ /WAV OPT, OBRA, SCALE ,RANGE,SQCN EXT, A,B, IDATE(30) , IHEAD(80) 
AA=21. 6*SQCN*RANGE*WAV** (- .3333333) 
EX-EXP (- EXT*RANGE) 
=(. 
RG2(1)=EX 
L=0 


4] 


or 


10 


15 
20 


DO 20 I=2,512 

IF (L) 10,5,15 

FF=976 .56* (I- 1) /SCALE 

ARG=AA* (FF**1 . 666667) 

IF(J.EQ.1) ARG= =ARGE (1. (FF*WAV/OPT) **. 3333333) 
IF(ARG.GT. 13 

IF((ARG-D) .L if L=-1 

RG? 1) =EX*EXP (- ARG) 

=ARG 


“ is 


Gta ose 


aT. 
RETURN 
END 


SUBROUTINE BITINV(K,L,M) 


C FINDS THE BIT INVERSION OF A (M) BIT NUMBER (J) TO (L) 
C THIS IS USED IN THE FOURIER TRANSFORM SUBROUTINE 
J=k 


— 


1 
2 
3 
4 


L=0 

DO 1 I=1,M 

L=L*2 

IF(J.NE.J/2*2) L=L+1 
J=3/2 


RETURN 
END 


SUBROUTINE CLOSE(IE1, IE2,1IE3) 

INTEGER*2 LU,LI,1H,IM,IS,IHS,IY,IMO,ID,SEL,VARY 
COMMON /O1/LU, LT, TH IM, 1S, THS, TY, INO, 1D 

CLOSE (LU, LOSTAT=IE1 

CLOSE(LI, IOSTAT=IE2 

CLOSE(9, IOSTAT=IE3) 

RETUR 

END 


SUBROUTINE WRITF (IER, BUF, L) 

DIMENSION BUF (512 

INTEGER*2 LU,LI,IH,IM,IS,1HS,1Y,IMO,1D,SEL,VARY 
COMMON /OI/LU,LI, IH, IM, IS, IHS, IY, IMO, 1D 

FORMAT (1H1) 

FORMAT(1X,/,1X,/,'''GAUSSIAN SOURCE FUNCTION! '') 
FORMAT (2X ,5E12.4) 

FORMAT (5X, 15, /,1X) 
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FORMAT(1X,/,1X,/,'''SOURCE LINE SPREAD FUNCTION''') 
FORMAT(1X,/,1X,/,'''FFT OF LINE SPREAD FUNCTION OF SOURCE''') 
FORMAT(1X,/,1X,/,'''OPTICS DIFFRACTION PATTERN 0(R)''') 
FORMAT(1X,/,1X,/,'''OPTIS LINE SPREAD FUNCTION OLS (Xr 
FORMAT(1X,/,1X,/,'''FFT OF OPTICS LINE SPREAD FUNCTION''' 
FORMAT(1X,/,1X,/,'''FFT OF SOURCE * FFT OF OPTICS''') 
FORMAT(1X,/,1X,/,'''LONG TERM MTF OF THE \TMOSPHERE' 

FORMAT 1K/, 1/41" FFT OF IMAGE LINE SPREAD FUNCTION LONG TERM''') 
FORMAT(1X,/,1X,/,'''IMAGE LINE SPREAD FUNCTION LONG TERM''') 
FORMAT(1X,/,1X,/,'''IMAGE POINT SPREAD FUNCTION LONG TERM''' 
FORMAT(1X/,1X,/,'''ENERGY INSIDE CIRCLE OF RADIUS R LONG TE 
FORMAT(1X,/,1X,/,'''SHORT TERM ATMOSPHERIC TPT) 
FORMAT(1X/1X/,'''FFT OF IMAGE LINE SPREAD FUNCTION SHORT TERM''') 
FORMAT(1X,/,1X,/,'''IMAGE LINE SPREAD FUNCTION SHORT TERM'!') 
FORMAT (1X,/,1X,/,'''IMAGE POINT SPREAD FUNCTION SHORT TERM! '' 
FORMAT 1X/1X/, "1 "ENERGY INSIDE CIRCLE OF RADIUS R SHORT TERM''') 
FORMAT (2X ,272,':',12,':',12,2X,12,'-',12,'-',14) 

IER=0 

WRITE(L 


J=(L+6 
GO TO (5579, elon tlm 23.2522 male sos Somo lias 
WRITE(LU, 2) 

GO TO 38° 

WRITE(LU, 10) 

GO TO 38 

neste) 

GO TO 


WRITE, 26) 
GO TO 
WRITE(LU, 34) 
GO TO 38 
WRITE(LU, 42) 
GO TO 38 
WRITE(LU,50) 
GO TO 38 
WRITE(LU,58) 
GO TO 38 
WRITE(LU, 66) 
0 T0 38 
WRITE(LU, 74) 
GO TO 38 
WRITE(LU,82) 
GO TO 


VRITE(LU, 90) 
GO TO 38 
WRITE(LU, 98) 
GO TO 38 
WRITE(LU, 106) 
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GO TO 38 
33 WRITE(LU, 114) 
GO TO 38 
35 WRITE(LU, 122) 
GO TO 38 
37 WRITE(LU, 130 
38 aoa 
RETUR 


20 WRITE LU,6) (BUF(T) ,I=1 64) 


6 FORMAT('WA = ',B12.5,' DIAMETER OF OPTICS = 
$E12.5,/,'OBSCURATION RATIO ' £12.5,' SCALE = 
$ £12.5,/, 'RANGE = ' F12.5,' CN SQUARED = 
$' £12.5,/, EXTINCTION COEEFICIENT = ' ,£12.5,' AMPLITUDE = 
$ ' F12.5,/,'SIGMA = '£12.5,/,2X,15A4, /,2X,20A4, / 
§ 2X, 20A4) 
RETURN 
END 


SUBROUTINE FFT(RG,ID,KA) 

C FOURIER TRANSFORM SUBROUTINE 
IMPLICIT COMPLEX (C) 
DIMENSION CE ae CG (1024) 
DIMENSION RG(512 
DATA IST ae .006135923 / 
CG(1) =RG(1 ) 

DOEIO. 1225512 
CG 1)=RG(1) 

10 CG(1026- 1)=CG(T) 
CG(513)=CG(512) 
IF(IST.NE.0) GO TO 2 
lito | 
DOS lai 
A=TH*FLOAT(I- 1) 
CE(T)-CERP(C 

CECORT) <0, 

Chi25 = One) 

DO 5 M=1,10 

aC iy 

L=2** (M- 1) 

DO 4 LL=1,L 

KC=LL- 1 

CALL BITINV (KC,N,9) 

CD=CE(N+1 

IF (ID. Eq. 1} CD=CONJG (CD) 

eZee kK 

KU=KL+K- t 

DO 3 T=KL,KU 

CA= C6(T)+CG( +k) *CD 

eer} =Ca =C¢6(1)- CG(1+k)*CD 
CA 


b> _— 


3 CG(I)= 
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CONTINUE 

CONTINUE 

DO 8 I=1,1024 

J=I-1 

CALL BITINV (J,N, 10) 
IF(E-N) 6, 6,7,8 

6 N=N+1 

CA= =06 (I) /(32. 0.) 


CG (I)=C6(N) 
CG(N)=CA 
CG (I)=CG(1) /(32. ,0.) 
CONTINUE 


DO 20 I=1,512 
RG (I)=CG ( (1), 
IF RoC) 'GT.1.E-6*RG(1)) KA=I 
20 CONTI 
RETURN 
END 


Ot 


CO =~] 


SUBROUTINE INT(K,TPOW) 
C CALCULATES THE FRACTION OF POWER INSIDE A CIRCLE OF RADIUS R 
C OR CALCULATES THE POWER INSIDE A CIRCLE OF RADIUS R 
COMMON /REAL/RG1 (512) ,RG2(512) 
COMMON / /WAV. OPT, OBRA, SCALE, RANGE, SQCN, EXT, A,B, IDATE(30) , IHEAD(80) 
an E*RANGE) **2 
RG2 1)= 7853982 12*RG2(1)*SR 
DO 3 1=2,512 
IF(I-K) 1,1,2 
RG2 (1) =RG2(I- 1)+6.28315E- 12*RG2(1) *SR*FLOAT(I- 1) 


RG2( (1) = =RG2(K) 
CONTINUE 

TPOW=RG2 (512) 

DO 5 I=1,512 
R62(T)= RG2(1) /TPOW 
RET 


END 


GOK be 


or 


SUBROUTINE LSF (RG3,K) 
C CONVERTS POINT SPREAD FUNCTION (R) TO LINE SPREAD FUNCTION (X) 
DIMENSION RG3(512) 
DO 2 I=1,K 
7= =R63(1) 


R-SORT (FLOAT((F- 1)**24J*J)) 
IR=IFIX(R 
RI=FLOAT (IR) 


RF=R- RI 
Z=Z+2*((1.-RF)*RG3(IR+1)+RF*RG3 (IR+2) ) 


— 
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KO 


10 
13 


el 
14 


12 
15 


20 
40 


1005 
60 


1140 
10 
1150 


J=J+1 
IF an ae K) GO TO 1 
RG3(I 

CONTI tr 

RETURN 

END 


SUBROUTINE OPEN 

INTEGER*2 LU,LI,IH,IM,IS,IHS,IY,IMO,ID,SEL,VARY 

COMMON /OI/LU,LI,IH,IM,IS,IHS,IY,1MO, 1D 

OPEN(7, IOSTAT=IE1,ERR=10,FILE='B: INPUT. DTA' ,STATUS='OLD' ) 
OPEN(8, IOSTAT=IE2,ERR=11,FILE='B: OUTPUT .DTA' ) 

OPEN 9° IOSTAT=IE3, ERR=12, FILE='B:0UT2.DTA') 

RETURN 

WRITE(*, 13) IE1 

FORMAT(' ERROR STOP IN OPENING INPUT FILE ERN = ',12) 
STOP 


WRITE(*, 14) IE2 


FORMAT(' ERROR STOP IN OPENING OUTPUT1 FILE ERN = ' 12) 
STOP 
WRITE(*, 15) IE3 

= '12) 


FORMAT(' ERROR STOP IN OPENING OUTPUT2 FILE ERN 
END 


SUBROUTINE OUTIN(BUF,J,K,L) 
INTEGER*2 NAME,LU,LI,IH,IM,IS,IHS,1Y,IMO,1ID,SEL,VARY 
COMMON /O1/LU,LL.1H, IM,15,1H8,1Y, INO, 1D 

DIMENSION BUF (1024) 

CALL GETTIM (1H,IM,1S,IHS) 

IF(J) 20,40,60 

CALL OPEN 

RETURN 

CALL READF(IER,BUF ,L) 

IF(IER.GE.0) RETURN 
WRITE(LU, 1005) IER, IN 

FORMAT ( ERROR STOP. READF ERROR. IER 


IF(K.LT.20)RETURN 
CALL WRITF (IER, BUF, L) 
rr lin 0) eet GO TO 70 


' 16, 'BLOCK # = ',I6) 


IF(IER.GE.0) RETURN 

WRITE(LU, 1140) IER, IN 

FORMAT ERROR STOP. WRITF ERROR. IER 
P 

WRITE(* ,1150) 

CALL OLSE (1 1, 1E2, IE3) 

FORMAT(' CALCULATION COMPLETE' ) 

STOP 

END 


' 16,' BLOCK # = ',I6) 
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SUBROUTINE ate BUF, iL) 

INTEGER*2 LU,LI,IH,IM,IS,IHS,IY,IMO,ID,SEL,VARY 
DIMENSION BUF (512) 

COMMON /OI/LU,LI,IH,IM,1IS,IHS,1Y,1MO,1D 
IF(L.EQ.1) GO TO 20 

ROMA F105) 

DO 10 T=1,512,8 


K=I+7 
ag ee 
NUE 


10 CONT 


) 
to © 


= WwW 


IER=1 
RETURN 

READ (LI ,2) (BUF (I) , 1=1,9) 
ee (Se, 5) 

READ(L 3) (B (I) ,1=10,24) 
FORMAT (1544 
READ (L 4) ( UF (1) ,1=25,64) 
FORMAT (2044) 

RETURN 

END 


47 


APPENDIX B 


A. MODULE INTERCONNECTIONS 


Hla 


Ce aan 


[Fig. 4] 
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ut 


WRITF 


WRITO 
CLOSE 


RESUL — 


Module interconnections 
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a ae | Gab feo 

NCO 
Yon) 

jem 

— 

<= 

(2 aa i 
ae Se aa 

i Se, 

oO Oo Ore 

Oo OO = 

> ae oe 

Ct | 

a = 
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APPENDIX C 
OUTPUT OF THE SAMPLE RUNS 


PeeOUTPUT OF RUN 1 


VARIABLE RADIUS AVERAGE USABLE POWER 

.8885E- 14 .240 .9770E+00 

.2030E- 13 .360  .4325E+00 

.3172E- 13 .470 .2527E+00 

.4313E- 13 .550 .1838E+00 

.5455E- 13 .630  .1395E+00 

.6597E- 13 700 .1125E+00 

. T738E- 13 770 =. 9253E- 01 

.8880E- 13 .830 .7925E-01 

.1002E- 12 .880 .7015E-01 

.1116F- 12 .940 .6115E-01 

.1231E-12 .990 .0483E- 01 

el S45E- 12 1.040 .4939F- 01 

.1459E- 12 1.080 .4553E- 01 

.1573E- 12 1.130 .4132E 01 

.1687E- 12 1.170 .3829E-01 

.1801E- 12 i210 © .3555E-01 

.1916E- 12 1.250 .3306E-01 

. 2030E- 12 1.290 .3080E-01 

.2144E- 12 LAS vAs, .2918E- 01 

.2258E- 12 1.360 .2726E- 01 

VARIABLE IS SQCN 

WAVELENGTH =  .10600E-05 DIAMETER OF OPTICS = .95000E-01 
OBSCURATION RATIO =  .00000E+00 SCALE =  .10000E+01 
RANGE =  .10000E+05 CN SQUARED = varied 
EXTINCTION COEEFICIENT = .50000E-04 AMPLITUDE =  .50000E+03 
SIGMA = .15000E+01 


05:18 19 12 1989 
Assumed Gaussian source for simulation 
Assumed optics information for run 1 


ol 


B.. OUTPUT OFSEUN 2 


VARTABLE RADIUS AVERAGE USABLE POWER 

. S000E- 05 .250 .1411E+01 

. 3105E- 04 .250 .1088K+01 

.OTILE 04 .250 .8383E+00 

.8316E- 04 .250 .6460E+00 

.1092E- 03 .250 .4978E+00 

.1353E- 03 290 =. 883 7E+00 

.1613E- 03 250 .2957E+00 

.1874E- 03 .250 .2279K+00 

.2134E- 03 .250 .1756E+00 

.2395E- 03 .250  .1353E+00 

.2655E- 03 .250 .1043E+00 

.2916E- 03 .250 .8037E-01 

.3176E- 03 .250 .6194E 01 

.3437E- 03 .250 .4773E- 01 

.3697E- 03 .250 .3678E-01 

. 3958E- 03 .250 .2835E- 01 

.4218E- 03 .290 .2185E-01 

.4479E- 03 .250 .1683E- 01 

.4739E- 03 290 .1297E-01 

. 5000E- 03 .250 .9998E- 02 

VARIABLE IS EXT 

WAVELENGTH =  .10600E-05 DIAMETER OF OPTICS = .95000E- 01 
OBSCURATION RATIO =  .00000E+00 SCALE =  .10000E+01 
RANGE =  .10000E+05 CN SQUARED = .10000E- 13 
EXTINCTION COEEFICIENT = varied AMPLITUDE = .50000E+03 
SIGMA = .15000E+01 


06:15 19 12 1989 
Assumed Gaussian source for simulation 
Assumed optics information for run 2 


52 





fee VUIPULT OF RUN 3 


VARIABLE RADIUS AVERAGE USABLE POWER 

.50008+04 .090 .2234E+01 

.5526E+04 105 .1952E+01 

.6053E+04 Pla .1901E+01 

.6579E+04 5 132 .1671E+01 

.1105E+04 .149 .1476E+01 

. 7632E+04 . 168 .1309E+01 

.8158E+04 .188 .1167E+01 

.8684E+04 . 200 .1136E+01 

.9211E+04 SPA .1016E+01 

.9737E+04 .243  .9120E+00 

.10268+05 .257 .8882E+00 

.1079E+05 281 . 7997E+00 

.1132E+05 .306 .7221E+00 

.11845+05 .3820 =. 703 3E+00 

.1237E+05 .346 .6368E+00 

.1289E+05 Bes .5781E+00 

.1342E+05 389 =. 56 30E+00 

.1395E+05 .418  .51245+00 

.1447E+05 .434  .4990E+00 

.1500EF+05 .465 .4551E+00 

VARIABLE IS RANGE 

WAVELENGTH = .10600E-05 DIAMETER OF OPTICS = .95000E-01 
OBSCURATION RATIO =  .00000E+00 SCALE = .10000E+01 
RANGE = varied CN SQUARED =  .10000EF- 13 
EXTINCTION COEEFICIENT = .50000E-04 AMPLITUDE =  .50000E+03 
SIGMA = .15000E+01 
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D. OUTPUT OF RUN 4 


VARIABLE RADIUS AVERAGE USABLE POWER 

.2500E- 01 .330  .5145E+00 

.3026E- 01 .310  .5840E+00 

. 3003E- 01 .290 =. 6672E+00 

-4079E- 01 .280 .7166E+00 

-4605E- 01 .280 .7173E+00 

.5132E- 01 .270 =. 7715E+00 

. 0658E- 01 .270 =. 7716E+00 

.6184E- 01 .260 .8321E+00 

.O711E-01 .260 = .8321E+00 

.1237E- 01 .260 .8321E+00 

. 1763E- 01 .260 =.8321E+00 

. 8289E- 01 .250  .9000E+00 

.8816E- 01 .250  .9000E+00 

.9342E- 01 .250  .9000E+00 

. 9868E- 01 .250  .9000E+00 

.1039E+00 .250  .9000E+00 

.1092E+00 .250  .9000E+00 

.1145E+00 .250  .9000E+00 

.1197E+00 .250  .9000E+00 

.1250E+00 .250 .9000E+00 

VARIABLE IS OPT(DIAMETER) 

WAVELENGTH =  .10600E-05 DIAMETER OF OPTICS = varied 
OBSCURATION RATIO =  .00000E+00 SCALE =  .10000E+01 
RANGE =  .10000E+05 CN SQUARED =  .100005- 13 
EXTINCTION COEEFICIENT = .50000E-04 AMPLITUDE =  .50000E+03 
SIGMA =  .15000E+01 


11:38 19 12 1989 
Assumed Gaussian source for simulation 
Assumed optics information for run 4 


o4 


jay 


GAUSSIAN SOURCE FUNCTION 


GRAPHICAL OUTPUT OF CALCULATION STEPS 


RE err ce 


500 


jam) & S & © —) © = 
wy © 2) om) Ww © ww 
a wv foe) oe) N N a — 


(1s/syem) AyIsuaiuy 


[Fig. 5] Gaussian source function 
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[Fig. 6] Source line spread function 
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(Fig. 7] FFT of source LSF 
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(Fig. 8] Optics diffraction pattern 
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[Fig. 9] Optics line spread function 
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(Fig. 10] FFT of optics LSF 
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(Fig. 11] FFT of source * FFT of optics 
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[Fig. 12] AMTF (long term) 
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FFT OF IMAGE LSF (LONG TERM) 
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[Fig. 13] FFT of image LSF (long term) 
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[Fig. 14] Image LSF (long term) 
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[Fig. 15] Image PSF (long term) 
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ENERGY INSIDE CIRCLE OF RADIUS R (LONG TERM) 
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(Fig. 16] Energy inside circle of radius r (long term) 
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(Fig. 17] AMTF (short term) 
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FFT OF IMAGE LSF (SHORT TERM) 
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(Fig. 18] FFT of image LSF (short term) 
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[Fig. 19] Image LSF (short term) 
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[Fig. 20] Image PSF (short term) 
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[Fig. 21] Energy inside circle of radius r (short term) 
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